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We note the similarity between BEC (Bose-Einstein Condensates) formed of atoms between which 
we have long-range attraction (and shorter-range repulsions) and the field theoretic " Q balls" . This 
allows us in particular to address the stability of various putative particle physics Q balls made 
of non-relativistic bosons (K 0, s, B°'s, and D 0, s) using variational methods of many-body NRS 
(Non-Relativistic Schrodinger) equation. 

PACS numbers: 



Introduction 

Phase transitions occur when a change of coupling or 
temperature makes the (free) energy of a new phase lower 
than that of the preexisting phase. In the field theo- 
retic formulation, this is manifest when the minimum of 
the effective potential U((f>), with the field <f> representing 
some order parameter, is shifted from cf> = to nonzero <f> 
value or to a degenerate manifold of such <j> values. In a 
simple mechanical analog the system represented by the 
single coordinate <j>, "rolls over" to the new stable min- 
imum. Some time ago, Sidney Coleman introduced[l] 
"Q balls" : classical field configurations stabilized by the 
global conserved charge (s) (Q) they carry. These corre- 
spond to configurations of, say, a charged field which are 
constant over a large region of space but vary in time: 
[<f> exp(— iwot)] and exp(nuot)] — corresponding to 

a constant charge density j° = p = (^0 + •(/) — + ^</>)/2i. 
The simple mechanical analog here is a system rotat- 
ing with uniform angular velocity in the <pi — (f>2 plane 
(with </> = <$>i + i 02). The "centrifugal" force gener- 
ated can then, under certain conditions specified below, 
make the representing "particle" come to equilibrium at 
a 4>q which is no longer a minimum of U ((/)). Coleman's 
suggestion of using the conserved baryon number as the 
global charge of the Q balls has been followed up in super- 
symmetric modcls[2]. Motivated by advances in BECs 
(Bose-Einstein condensates) , we investigate in this paper 
their relation to Q balls and use variational NRS (non- 
relativistic Schrodinger) equation methods to prove the 
stability against strong interaction decays of Strangeness, 
Charm, and Beauty balls. 

Field theoretic Q balls are more general than the non- 
relativistic limit on which we focus here. Thus [2] the 
above-mentioned baryonic Q balls, made of squark con- 
densates, are stable against "weak interaction like" de- 
cays of the heavy individual squarks by having the Q balls 
very tightly bound with masses proportional to a frac- 
tional power of the total baryon number iV*- 1-6 -* rather 



than N 1 as expected for non-relativistic weakly bound 
matter. 

Still, the equivalence — in some limits — of the field- 
theoretic and many-body descriptions of the same "con- 
densations" is interesting and helpful. 

There are two distinct, though interrelated, types of 
condensation. BECs obtain when a large number, N, of 
(bosonic) atoms are trapped within a radius R and then 
cooled down to nano-Kelvin temperatures. The BE con- 
densation is in momentum space: a finite fraction of the 
atoms are in the lowest mode of the trap corresponding 
to p = 0. This manifest when the trap is suddenly re- 
moved by having these p = v = atoms hardly move. [3] 
Boson-boson interactions modify the BEC but are not re- 
sponsible for BEC phenomenon in the first place. BEC is 
best understood when the atoms are non-interacting[4]. 
On the other hand the coordinate space "condensation" 
of bosons into Coleman's "Q balls" is due to attraction 
between the bosons. Thus in the field theoretical for- 
mulation Coleman has proved that stable Q balls exist 
if and only if the "potential" U(<fi) in the effective low- 
energy Lagrangian for the system satisfies: 

(i) — has a minimum lower than — (1) 

q) 2 - 2 

with /i the mass in the "free" part of U(<j>) : ^(f) 2 ■ 

The size R of the spherical Q ball and its density n = 
(3Q)/ (4ttR 3 ) are fixed by the overall Q and (f>o — the field 
for which the above minimum is achieved. 

Condition (i) implies over all attractive interactions 
between the 4> bosons, say, the kaons in strangeness balls, 
at an appropriate density. In particular, (i) holds if the 
coefficient of the lowest (</>) 4 term in U{4>) is negative: 

„2 

(ii) U(4>) = —(f) 2 — \(j) 4 + higher order terms 

(2) 

corresponding to attractive S-wave scattering length (A > 
is implicit). Note that U{<j)) is an effective Lagrangian, 
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which is not used in loops inside Feynman diagrams. 
Therefore U (0) can have (and indeed has) higher-order 
positive non-renormalizable (</> + 0)" terms-ensuring a fi- 
nite </>o and a spectrum which is bounded from below. 

While (ii) — > (i), namely, an attractive S-wave scat- 
tering length (negative </> 4 coefficient), implies stable Q 
balls, (ii) is not required. Thus, a nontrivial minimum, 
<po, of U (<fi) / (4>) 2 can obtain with a negative (0 + </>) 3 term 
overcoming at some <fi the positive ^{fp) 2 + -M0) 4 . Still 
this can be problematic: the minimum of £/(0)/(</>) 2 may 
now be at a large 4>q, say, </>o >> fj. where the NR many- 
body approach that we want to compare with next, fails. 

Note that Condition (ii) alone, without knowledge of 
the higher-order terms, does not fix the size (i?) or the 
density (n) of the Q balls: with only the [— A0 4 ] term 
present both cf>o and n are infinite! 

Forming a Spatial Droplet of Non-relativistic 
Bosons: The Many-Body Approach 

Let N non-relativistic identical bosons of mass m in- 
teract via potentials V{\r,i — Tj\). The basic question we 
address is: "For which potentials the NR bosons coalesce 
into "Q balls" with nonzero density when N — > oo ?" 

For finite range potentials the N-body Schrodinger 
equation if I*) = E\$) with 

^ = Ef|+E[^-^)] (a) 

is trivially solved by \I r which is a product of single parti- 
cle wave functions, each of which is (approaching) a con- 
stant. The spread-out particles have vanishing kinetic en- 
ergies and vanishing mutual interactions and thus E = 0. 
For many potentials a lower (negative) energy state ex- 
ists with the N bosons in a sphere of radius R. Using 
the basic variational principle we derive next several suf- 
ficient conditions for that to happen. In Sees. Ill, IV 
below we argue that some of these sufficient conditions 
are met in the case of many K°'s, D 0, s and B 0, s which 
therefore will make droplets or Strangeness, Charm and 
Beauty balls which are strong interaction stable. 

Let us first use a simple trial wave function with all 
bosons in the same state — the ground state of a large 
radius R spherical cavity: 

** = IJWVO)- (4) 

i 

The expectation values of the kinetic and potential ener- 
gies appearing in (^tlH^t) = (K) + (V) are: 

(K) ~ Nh 2 /(2mR 2 ) (a) (5) 

and, if R >> tq = the range of the potential, (V) ~ 
[Nnv] - [N 2 /(4irR 3 /3)}v, with 

v = 4tt / V(r) r 2 dr {b) (5) 



and n=N/ (volume) the particle number density. A pos- 
sible relation to the previous section stems from the fact 
that, up to kinematic factors, v is the Born approxima- 
tion for the S-wave boson-boson scattering length. For 
dilute systems with inter-particle separation which far 
exceed the range of the potential: 

d = n- ( V 3) »r , (6) 

only the integrated potential v effects the threshold scat- 
tering. The threshold scattering amplitude is therefore 
reproduced also by a "pseudo-potential" of a local delta 
function form: vS 3 (r). Such a potential roughly cor- 
responds to the negative local ~\(<f> A ) term — which in 
the field theoretic formulation suffices to ensure stable Q 
balls. To complete the analogy with the NRS case we 
show that if v < 0, then also the expectation value of the 
energy (if) < 0. This readily follows from the different 
scalings of the (expectation values of) the positive kinetic 
and negative potential energy N, the number of bosons 
in the system: (K) ~ NW « ~ N for TV -> oo. 
By the variational principle the energy of the true N- 
body ground state is lower than (if) and also negative 
and a many-body bound droplet or Q ball stable state 
indeed exists if: 

(tfJ.^/vW^^O. a "droplet" state exists. 

(7) 

Like (ii), (ii') does not fix the actual size/density of the 
"droplet". Still the two conditions are not equivalent. 

The quantity v, depends only on the potential V(r) 
and not the mass m. It is (proportional to) the actual 
S-wave scattering length (or to the coefficient of the </> 4 
term in U (4>) only in the Born (or dilute system) approx- 
imation. 

Scattering theory [5] implies that the S-wave scatter- 
ing length is attractive if and only if the NRS (non- 
relativistic Schrodinger) two-body system has bound 
states. Thus the NRS equivalent of the field theoretic 
condition (ii) is having a two-body bound state of the 
NRS equation. Indeed as we directly show below hav- 
ing two-boson bound states guarantees N bosons bound 
state. Note, however, the independence of the two 
NRS criteria: The criterion (ii') ensuring N — > oo NRS 
"droplets" does not ensure an S-wave two-body bound 
state. In three dimensions the latter requires not only an 
"attractive" potential, but also sufficiently strong attrac- 
tion. 

Various criteria for V(r) to have bound states in a 
NR two-body system with reduced mass m exist. [6] Yet 
there is no general if and only if criterion, short of solving 
the Schrodinger equation. Finding if a bound state of 
N — > oo bosons exists need not be easier. Indeed the field 
theoretic Criterion (i) requires the full effective potential 
U{4>). The coefficients in the power series for the latter 
are the threshold scattering amplitudes for any number 
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of particles and cannot be found short of solving exactly 
the field theory. 



"Strangeness-Beauty Balls" and the 
Non-Relativistic Schrodinger Equation 

For many atomic and other systems, the many-body 
NRS treatment preceded field theoretic approaches and 
Laughlin's celebrated variational wave function for the 
fractional quantum Hall effect is a prime example of this. 
In recent years much effort has been devoted to applying 
field theory in general, and effective field theories and 
effective Lagrangians in particular, to such problems and 
this has been also the case for the Q balls. [7] 

Chiral Lagrangians were used to check if Coleman's 
Criterion (i) holds for the K° system. These Lagrangians 
involve higher derivative terms and are fixed by an over- 
all fit to data yet they do not offer much intuitive un- 
derstanding. Here we follow the reverse program of "de- 
mystifying" Q balls and trying to explain them — at least 
in the NRS regime — as simple spatial condensations of 
many bosons with appropriate potentials. Strangeness- 
balls, and even more so, Charm/Beauty-balls with den- 
sities n < m K 3 are NR. In all cases we expect to have 
essentially the same boson-boson potentials as the later 
are controlled by the common light d quark. Hence, us- 
ing K° - K° potentials, V K (r), etc., to find if U K°, etc., 
Droplets" form — which we do next — is justified. 

A first key observation is that the K — K potential, as 
that between any two identical, neutral, pscudo-scalars, 
is attractive at "large" ~ Fermi distance. To show this 
we write V(r) as a superposition (integral) of Yukawa 
potentials due to all exchanges [8]: 



V(r) 



-fir 



(8) 



Parity conservation forbids a KKix vertex and the 
lightest exchanged system controlling V(r) at r — > oo 
is that of two pions. Further, the t = component dom- 
inates at the 7T7T threshold: 

a{p 2 ) near fx = 2m{n) ~ \f(l = 0)(K + K —t 7T7t)| 2 (9) 

With f(£ = 0), the tttt <-> KK S-wave amplitude. This 
expression is clearly positive and recalling the minus sign 
in the definition of V(r) we find that the longest range 
two-pion exchange potential is indeed attractive. This 
can be also directly shown to be the case for the exchange 
of a ++ state which can be an S-wave resonance in the 
7T7T S-wave system. The above reasoning is similar to that 
used in [8] to derive the well-known attractive Casimir 
Polder (i.e retarded van der Waals) and the regular van 
der Waals two photon-exchange potential between two 
identical neutral atoms. 

At short distances of the order of the size of the sd com- 
posite state, which is the kaon in quark modcls/QCD, the 



KK potential becomes repulsive. As in atomic physics 
this is due to the Pauli principle — operating here between 
the identical d or s quarks in the two K 0, s. The repul- 
sion can also be viewed as being due to the exchange of 
the u) (p) vector mesons: The two K°'s have the same 
hypercharge (isospin) to which u> (p) couple [9]; no light 
++ , scalar "er" meson has been established. Yet, the 
exchange of such an entity with appropriate mass and 
coupling to nucleons g(o~; NN) (which may represent the 
box diagram with two pions exchanged and intermediate 
N and A(1230) states) could, along with u exchange, 
dominate nuclear binding [10]. 
Parameterizing, 



V K (r) = -g(af 



-(m a r) 



9{Vf 



(10) 



with V = ui or p at a common mass m(V) and g(V) 2 the 
sum of the (squared) p and lo couplings to kaons we find: 



v K 



m(o~) 5 



m(vy 



> 



(11) 



To evaluate (II) we take g(p,KK) = g(us,KK) ~ 
(l/S)g(w,NN) and g{a,KK) ~ (l/3)g(a,NN), as 
suggested by counting the numbers of non-strange 
quarks. Using the values of g(uj 7 NN) 2 /rriy and 
g(cr, NN) 2 /m(a) 2 suggested by fitting nuclear matter 
in Eq. (14.27) in [10] we find that v K ~ (-2.35 + 
3.45)/m^r > 0, so that Criterion (if) is not satisfied. 
The box diagrams with K* intermediates for KK sug- 
gests a g 2 KK which is somewhat bigger than the previ- 
ous [(l/9)<?(er, NN) 2 ] estimate. Also cutting off V K (r) of 
Eq.(10) in evaluating vk at, say, 0.3 Fermi — suggested 
by the fact that the K° is composite at such a scale and 
can no longer be treated as a point source of the a, u> 
fields — further reduces the repulsive relative to the at- 
tractive contribution and a negative vk is not excluded. 

Still vk > is likely and we face the question: "Does a 
positive vk exclude stable K° droplets?" This is the case 
if we insist on N-body wave functions which are products 
of N identical one-particle wave functions ip{ r i)- How- 
ever, including the (Jastrow) product of N(N — 1)/2 two- 
particle functions 

*(trteJ) = IJW r *)) II (/(ki-^-l)) (12) 

• (*>i) 

we can have an N-body bound state even if v > 0. 

To illustrate this we fix the potential Vff(r) (and vk) 
and increase the mass. (Taking m = ms ~ Wvtik corre- 
sponds, in the approximation of universal potentials be- 
tween pairs of heavy mesons — and treating s as a heavy 
quark, to discussing "Beauty" (rather than Strangeness) 
balls. While vk ~ vb ~ v the system can now have even 
two-body bound states. 

In the m — ► oc limit these are localized around the min- 
imum of V(r) at r = ?*o, forming a "Vibrational Band" 
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with spacings ~ ([^"Ml 7 * = J'oJ/m)^ 1 / 2 ^ . Once the po- 
tential V(r) has two-body bound states, N-body bound 
states are guaranteed. This is verified by using in Eq. 
(12) f(r) = i' {r) with r = \n - r 3 \ satisfying: 

(- — TT + V(T))Mr) = EoMr) (13) 

with Eq < the negative energy of the two-body bound 
state, and operating on the above trial function with the 
N-body Hamiltonian of Eq. (3) above. (Here, the re- 
duced pair mass m* = m/2 replaces the single particle 
mass in the Schrodinger equation.) Even when the po- 
tential is too weak relative to the kinetic term to have a 
two-body bound state, using Eq. (12) with f(r) of the 
above general form, namely, peaking at the minimum of 
the potential and being small at small r's where V(r) is 
repulsive, lowers (H) relative to its value for the product 
of one-particle functions. As we show in some detail in 
the next section, this can yield the desired N-body bound 
state even when Condition (if) fails, and also there are 
no two-body bound states. In this case we have a fully 
symmetric N-boson bound state with the f{\fi — rj\) fac- 
tors peaking at r = r$, and determining the density of 
the N-boson droplet to be: 

n~([(47r)/3]rg)C- 1 ) (a) (14) 

or, equivalently, the radius of the droplet 

R^NWiro (b) (14) 

This is reminiscent of U(<f>) / (4>) 2 and its nontrivial min- 
imum at 4>o — fixing the radius and density of the Q balls 
in the field theoretic formulation. There the nontrivial 
cf>o obtains via the interplay between a negative <fi 4 term 
and positive higher-order terms. In the present NRS case 
the minimum at r$ reflects attraction (repulsion) at long 
(short) ranges. The higher (</>)" terms are prominent 
at large densities — just like the strong short-range repul- 
sions in NRS. 

Still, a U{4>) <-> V(l/r) analogy is rather limited: The 
effective potential U derives from the fundamental La- 
grangian of the field theory, say, QCD for the above cases, 
whereas the potential V(r) is the primary entity in the 
NRS approach. A closer analog of U{<j>) is the derived 
quantity [E/N](n) — the energy per particle for a given 
density[ll] in NRS. E = E[N; R] is the ground state en- 
ergy of the Hamiltonian in Eq.(3) of N = h[(4tt)/3] • R 3 
bosons uniformly distributed (after averaging over corre- 
lations) in a sphere of radius R. When N and R tend 
to oo keeping n fixed, a stable droplet of density n ob- 
tains that if and only if the minimum of [E/N](n) is at 
< no < oo and is negative [12]. 



Binding and BECs in the Presence of Strong, 
Short-Range Repulsions 

Short-range repulsions do not hinder BEC for dilute 
atomic systems in traps. 

Let N-bosons be in the trap and add one more. To 
see the issue most clearly, assume first that the N-bosons 
are "frozen" at specific locations r° inside the trap. The 
added k « N bosons will be in the ground state of the 
total potential: 

V = V (trap) (r)+ V(\ r i-r\) (15) 

(i=l,..iV) 

— provided that the potential 15 can bind a particle of 
mass m (which is clearly the case for an attractive V( trap ) 
and sufficiently dilute atoms). 

Conceivably such a setup, of interest in its own right, 
can be experimentally realized. Let the trapped N atoms 
form a 3-D lattice generated by a standing wave pattern 
of three lasers. Let the added atoms be of a different 
species which interact with the first TV atoms via the 
potential V in Eq. (15) above. In particular, we need to 
choose a species which is almost unaffected by the laser 
fields. 

"Freezing" N out of N + 1 is artificial and adding one 
extra boson causes each of the previous TV-bosons to ad- 
just by order 1/N, modifying the binding energy by 0(1). 
A key point is that the adjustments are likely to lower 
the energy and neglecting those is appropriate if we only 
want to verify that the extra particle binds. 

As the system becomes denser, stronger two-boson cor- 
relations and higher momentum components build up. 
The completely symmetric N-body state, while no longer 
factoring into independent N single-particle wave func- 
tions, still exhibits coherence features unique to BECs. 

Our main interest however is not this, but rather the 
"droplet formation" problem posed in the previous sec- 
tion. To address this problem we use Eq. (15) without 
V(trap) an d two-body potentials which are attractive (re- 
pulsive) at long (short) distances (Fig 1.) 

A likely configuration the N particles r° is at the 
vertices of a simple cubic lattice of lattice constant 
d = n^ 1 / 3 '. A simplified version of the problem which 
serves as a criterion for formation of a droplet of density 
n = d~ 3 is: 

(iii) "Does a particle of mass m bind to a cu- 
bic lattice with lattice constant d and common po- 
tentials V(\r — r.i\) centered at all lattice points 
r .?» 

— a problem which may also be useful in discussing 
trapping of light in some "dielectric lattices" via the 
Helmholtz equation. 

It is difficult to address it in the most general case, yet 
the following suggests that bindings are likely even when 
a single potential V(r) fails to bind and/or to satisfy v < 



FIG. 1: The KK, DD or BB NRS potential of interest. It 
is repulsive at short distances r < n and attractive for r > 
ri having an asymptotic - — -g — behavior due to two-pion 
exchange (qualitatively similar to a Lennard- Jones potential). 
The minimum of the potential is indicated by the dotted line 
is at r . 

0. Let ro, ri with ro > r\ be the points indicated in Fig. 1 
where the common radial potential V(r) has its minimum 
and where it changes sign, being repulsive (attractive) 
for r < ri (r > ri). Consider a unit lattice cell of side 
of length d with eight potential centers at its corners. 
We focus first on the "dilute" case with d » 2r > 
2r\. Apart from the eight spherical octants of radius 
r\ at the corners, V of Eq. (15) above is attractive at 
all the remaining part of the unit cell. The longer-range 
tails of the other potential centers make for an attractive, 
negative contribution in the form of a "Madelung sum" 
at any, (say, corner) point. This increases the volume of 
the connected region within the unit cell where V (sum) 
< beyond the minimal value: 

4-7T 

d 3 - —ri [~ -93 d 3 if d > 4r*i !] . (16) 

Thus we find that the particle is inside a ( very loose!) 
cage of size ~ d 3 "cornered in" by the repulsive potentials 
centered at the eight corners of the unit cell of the cu- 
bic lattice — as in a three-dimensional analog of a carton 
egg holder. Actually, in the dilute case considered here, 
the particle can "roam" over all the lattice, going into 
neighboring cells via the large circular openings of radius 
d/2 — ?*i between neighboring cells, lowering the energy. 
To show this more clearly we simplify the problem by: (a) 
replacing the repulsive potentials within spheres of radius 
7*1 centered at the vertices of the lattice by infinite, pos- 
itive, "square-well" potentials within the circumscribing 
cubes of size 2ri centered — just like the above spheres — 
at the corners of the lattice cell, and (b) replacing the 
attractive potential within the remaining region by their 
volume average, 

-uo = v-/d 3 (17) 



FIG. 2: A schematic picture of the attractive and repulsive 
regions (after step b has been implemented) inside the unit 
cell — the overall cube of side d. The repulsive regions (with 
V = +oo) are the eight smaller cubes of side 7*1 at the eight 
corners of the unit cell — the big overall cube. Three of these 
are shown as the full line small cubes in the figure. The po- 
tential is attractive (negative) outside these eight cubes with 
an averaged constant value. Using additional dashed lines we 
also illustrate two out of the twelve rectangular parallelepiped 
(RP) along the edges of the cube, between a pair of adjacent 
repulsive boxes. These can be viewed as "wave guides" for the 
particle in the potential where we have vanishing boundary 
conditions on the small (rf in area) sides of the RP and the 
particle is propagating in directions perpendicular to the long 
(d — 2 n in length) sides of the RP eg., in the top-front RP the 
wave guides in the up-down and forward-backward directions, 
and in the right-most RP illustrated in the up-down and the 
left-right directions. 

tential. (Sec Fig. 2.) 

Step (a) increases the volume of the region where the 
potential is repulsive by 6/ir and the value of the poten- 
tial which vanishes at its boundary to 00. This clearly 
increases the energy of the particle in the lattice. It is less 
obvious but still true that performing step (b), namely, 
replacing the attractive part of the potential by its av- 
erage value, also increases the energy of the particle in 
the lattice. To see this, note that the true ground state 
wave function for the original potential tends to con- 
centrate away from the repulsive regions and thus can 
sense better the attractive potential prevailing around 
the center of the above "cage", i.e., unit lattice cell, and 
this attractive potential is stronger before the averaging 
is performed. Hence, if after performing steps (a) and 
(b) above the particle still binds to the lattice, namely, 
(E (ground) < 0), then even deeper bound states are likely 
for the original problem. 

To test for a bound state we estimate the expecta- 
tion value of the Hamiltonian (H) = (V) + (T) in the 
state (ground))- Since the (normalized) wave function 
is nonzero only outside the repulsive regions, the expec- 
tation value of the constant attractive potential is: 



with v , the integral over the attractive part of the po- 



(V) = -u = v~/d 3 . (18) 
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The kinetic energy stems from the boundary condi- 
tions: tp = at the boundaries of the repulsive cubes. 
The twelve regions, of volume r\{d—2 r±) each, inside the 
d 3 unit cell between pairs of repulsive corner cubes can be 
viewed as short section parts of "wave guides" . The ki- 
netic energy can be roughly approximated as that of the 
lowest mode of the "waveguides": T = h 2 /[2m(d— 2ri) 2 ] 
and its expectation (T) is weighted by the fraction of the 
unit cell, / = \2r\ (d — 2ri)/d 3 , occupied by the "wave 
guide" sections, namely 



(T) = }T 



6r 2 



m d 3 (d-2n)' 



(19) 



The condition for binding, (T) + (V) < then becomes 

h 2 6r 2 



> 



m d — 2r\ 



(20) 



Since v is independent of d, the last condition can 
clearly be satisfied for d >> r\. 

As we gradually decrease d (relative to the distance 
scales ro and T\ of the potential V(r)), the average value 
of the attractive part, namely, \uq\ of Eq. (18) above 
increases, thereby enhancing the binding. 

However, once d ~ r± the "cages" trapping the particle 
within each unit cell become tighter and the particle can 
only tunnel between the different cells. The energy then 
rises and the bound state disappears. 

Finding the optimal d (or density of the droplet n) 
and the corresponding binding: [E/N](n) at this num- 
ber density requires detailed calculations beyond what 
we have attempted here [13]. 

Showing that one extra particle can be bound in a 
periodic box of size L where the previous TV are located at 
prescribed positions r ® (as we did above when the r ° were 
the nodes of a regular simple cubic crystal) is only one 
step towards proving the existence of large Kaon droplets 
with the specific KK potential above. Indeed: 

i) While we optimize the inter-particle separation d ~ 
7V( -1 / 3 ) • L to minimize ejy, the energy of the TV + 1 th 
particle added to the lattice, we should verify that the 
same d also allows each of the original TV particles to 
bind around the empty lattice site that it occupies. 

ii) We should show that the particular frozen r° ar- 
rangement at all lattice sites within the periodic L 3 box 
represents the "worst case" for binding the TV + 1 th par- 
ticle and any rearrangement of the TV bosons inside the 
box allows stronger binding of the TV + 1 th boson. 

iii) Using i) and ii) above we then show that enlarging 
the system from TV to TV + 1 bosons and simultaneously 
letting the cube size adjust to the new optimal length, 
L(N) — ■> L(N + 1), the energy is lowered by more than 
the above binding: 

E[(N + 1), L(N + 1)] < E[N, L(N)} + e N (21) 



This last step is readily achieved by employing 
E[(N + 1),L(N +1)} = 
Y[d 3 r° [ d 3 r N+1 V(r N+1 ;r°) 



x[( E 

L—1....N 

h 2 



2m 



V 2 + £ V (rlr°)) 



— V 2 



i>j = l...N 

£ V(rlr (N+1) )] 



-I...N 



*(nv+i;r ? °). 



(22) 



We use the suffixes on the first N coordinates to 
emphasize that in evaluating the expectation of the ki- 
netic energy and N interactions of the N + 1 th parti- 
cle by doing the innermost r^+i) integration, these first 
coordinates are "frozen". This integral then yields 
e[r°, .., rjy-], the binding energy of the extra particle to 
the first frozen TV which is then further averaged over all 

using the normalized measure provided by the density 
function of the first TV particles. If (ii) holds this yields 
a value smaller than than ejv above corresponding to the 
special case of a perfect full lattice of size L 3 (N + 1). 

Next consider the first term in the last integral, namely, 
H(N), the part of the Hamiltonian pertaining to the first 
TV bosons. Since H(N) does not depend on rny+i), the 
Normalized integral over H(N) simply yields E[N, L(N+ 
1)] which exceeds E[N, L(N)]. Adding the two terms we 
find at the desired inequality (21). 

Summing these over n < TV we yield a lower bound on 
the binding of the droplet; 



E N < Ne 



N 



(23) 



For a potential which is an attractive constant apart 
from hard core cubes of size 2r\ around each boson, i.e., 
the case studied above this last inequality, is directly 
proven in one dimension in Appendix I. Unfortunately, 
this particular elegant method does not readily generalize 
to three dimensions. 

Coming back to issues i) and ii) , clearly filling up com- 
pletely a simple cubic lattice with the previous TV bosons- 
rather than leaving one "hole" i.e., a vacant lattice site 
where the extra particle nicely fits lowers the volume of 
the "cage" in which the TVth particle is free to roam. 
This, in turn, elevates its kinetic energy. The attractive 
potential energy is {8V[(V3/2)<i] + 24U[(yiT/2)d] + ...}. 

Compare this to the "optimal" regular arrangement 
with one vacancy in the regular cubic lattice. The free 
volume where no strong repulsion occurs is now ~ twice 
as large making roughly for 2~( 2 / 3 ) ~ .65 times lower 
kinetic energy, whereas the attractive potential here is 
6V[d] + l2V[dV2] + .... For d ~ r 2 - 2ino~ ^ the latter 
is similar to the potential energy in the previous case. 
Hence we find that certain d values which allow binding 
of the N + 1 th particle to a perfect lattice make for even 
stronger binding inside the lattice. 
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Next let us consider a random rearrangement of the r° 
inside the L 3 cube with the same average number den- 
sity Rd~ 3 . It will contain pairs, triplets etc., of parti- 
cles which are nearer to each other than d and this will 
be compensated by having nearby "lacunas" of under- 
dense points. The latter constitute ideal placements of 
the extra TV + 1 th particle: It will have more free space 
to move and at the same time will be more strongly at- 
tracted to the dense "clusters". Note that in evaluating 
e[r°, rj, r^f] the binding of the extra N+l fh particle to 
the frozen TV particles at r?, we need not worry about the 
fact that the "crowding" of some of the frozen vertices 
raises their mutual interaction energies . 

Some Concluding Remarks 

In this paper we have shown a close correspondence be- 
tween the field theoretic concept of Q balls and droplets 
of coherent NR bosonic matter and elaborated at some 
length on the criteria, in a NRS picture, for forming the 
such droplets. 

While recently field theoretic/effective Lagrangian 
methods — such as those used by Coleman in predicting Q 
balls — are broadly applied to many-body physics, we find 
that the traditional variational NRS approaches can be 
used to prove the existence of some "particle physics" Q 
balls. Unfortunately these Q balls with Q = Strangeness, 
Charm and Beauty, while stable against decays via strong 
interactions do decay rather quickly via weak interac- 
tions. A single K° decays in ~ 10~ 10 sec (and D°,B° 
decay 100 times faster). A droplet of TV non-relativistic, 
weakly bound neutral bosons will start disintegrating af- 
ter times which are 1/TV shorter than the decay time of 
a single boson. If the minimal number of K 0, s can be 
viewed as an TV — > oo "droplet" is ~ 100, we need to 
assemble within ~ (5 Fermi) 3 100 slow K°'s in a picosec- 
ond!, which seems impractical. 

We would like to note, however, that the general fea- 
tures of the KK potential facilitating droplet formation, 
namely, attraction at "long" ~ (Fermi) distances and 
repulsion at "short" 1/3 Fermi distances, hold also for the 
if-Nucleon (and the TV — TV!) system. While there are 
both few body and "Droplets" of nucleons (a.k.a. "Sta- 
ble Nuclei" ) no K — TV bound states exist and the K — TV 
S-wave scattering length is known to be repulsive. 

However K° droplets exist despite the absence of a 
KK bound state. Could a K° which is free to move in a 
pre-existing large nucleus (so long as it avoids getting too 
close to the nucleons) bind to the latter? This is clearly 
not evident since the nuclear density and size are fixed by 
nucleon-nucleon interaction and not by K — TV physics. 

It is amusing to note, however, that the existence of 
such states would manifest experimentally as follows: Let 
a beam of slow K + charge exchange on a heavy nucleus 
with the resulting K° binding to the nucleus. The sub- 



sequent K a — > 7r + 7r~ decay ~ 10~ 10 sec later is likely to 
break the nucleus, yielding a spectacular "Star" which, 
in the absence of a A'°-nucleus bound state, should not 
occur. 

APPENDIX I 

In the "frozen TV" variant, the TV + I th particle is re- 
stricted to the TV intervals of size L' /TV between pairs 
of existing particles (periodic boundary conditions avoid 
half size end intervals) with L' = L — 2Nr\ the effective 
length allowed by the constraints. There is no tunneling 
between these TV intervals and the energy is : 

e N = H 2 /[2m(L'/N) 2 } (24) 

Note that Ne^ is the energy of one fictitious represen- 
tative particle moving in an TV dimensional cube of 
side L' /TV. The no-tunneling rigidity which is an arti- 
fact of one dimension, reflects also in the full TV-body 
problem of finding the ground state of the Hamiltonian 
S ~§^z f° r N particles in the (0, L') interval via the 
impenetrability — or ordering condition — 

< xi < x 2 < —-Xn < L' . (25) 

Equivalcntly we have one particle restricted to the above 
TV dimensional parallelepiped (Pn) satisfying the free 
Schrodinger equation in TV dimensions. 

(The TV! parallelepipeds P/v's obtained by permuting 
the Xi have the same volume L' N /Nl and their union 
constitutes a cube of side L'). contains the TV- 

dimensional cube Cn of side L/TV : 

< x 1 < L/N, 
L/N <x 2 < 2L/TV, 

(jfe - 1)L/N <x k < kL/N, 

(TV — 1)L/N <x N <L, 

and IIPjvII > ||Cjv||. (26) 

The variational principle and the above containment re- 
lation imply a lower energy for the full problem. 
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